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We compute the energy spectrum of the ground state of a 2D Dirac electron in the 
presence of a Coulomb potential and a constant magnetic field perpendicular to the 
plane where the the electron is confined. With the help of a mixed-basis variational 
method we compute the wave function and the energy level and show how it depends 
on the magnetic field strength. We compare the results with those obtained numerically 
as well as in the non rclativistic limit. 



1. Introduction 

Over the past years the study of systems of non-relativistic electrons confined to 
a plane in an electromagnetic field background has attracted much attention in 
view of possible applications. This problem is of practical interest because of the 
technological advances in nanofabrication technology that have made the creation 
of low dimensional structures like quantum wells, quantum wires and quantum dots 
possible. 1 The relativistic extension of this problem has also turned out to be of 
importance in the description of quantum two-dimensional phenomena such as the 
quantum Hall effect and high-temperature superconductivity. 2 Different condensed 
matter physics phenomena point to the existence of (2+l)-dimensional with an en- 
ergy spectrum determined by the Dirac equation Hamiltonian. 3, In particular, the 
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degenerate planar semiconductor with low-energy electron dynamics is assumed to 
admit an adequate description in terms of the (2+l)-dimensional relativistic Dirac 
theory. 5 In conclusion, the study of physical effects occurring in (2+1) dimensional 
systems of charged particles in strong external fields is an interesting problem from 
the theoretical point of view as well from its practical applications 6,7 ' 8,9 ' 10 . 

In order to analyze relativistic quantum effects in the presence of strong electro- 
magnetic fields one should be able to compute the Green function or to find exact 
solutions of the Dirac equation. Regrettably, the Dirac equation is exact solvable 
only in a very restricted family of electromagnetic configurations. 11 

The non relativistic two-dimensional Hamiltonian describing the Coulomb in- 
teraction — jr, between a conduction electron and donor impurity center when a 
constant magnetic B field is applied perpendicular to the plane of motion has been 
thoroughly discussed in the literature. 14,15,16,17,18 Despite its simple form their so- 
lutions can not be expressed in terms of special functions. An analogous situation 
occurs when we deal with the (2+1) Dirac equation, therefore one has to apply 
numerical and approximate techniques in order to compute the energy spectrum 
and the corresponding wave functions. 

In order to compute the energy spectrum of the relativistic hydrogen atom we 
have at our disposal different numerical and semi-analytic methods. The Kato and 
and Lehman-Maelhy methods have also been applied to solve the Dirac equation 
without variational collapse 12 A large finite-basis-set method has been applied in 
order to calculate in 3D relativistic and nonrelativistic binding energies of an elec- 
tron in a Coulomb field a magnetic field. 13 

It is the purpose of the present article to compute the energy spectrum of the 
ground state of the (2+1) Dirac equation in the presence of a Coulomb interaction 
and a constant magnetic field perpendicular to the plane where the electron moves. 
For this purpose we use a two-term mixed-basis variational approach. 16 We also 
compute the non relativistic limit with the help of the Pauli equation. Finally, we 
compare the results with those obtained numerically with the help of the Schwartz 
method, 19 which is a generalization of the mesh point technique. 

The article is arranged as follows: In Sec. 2 we write the (2+1) Dirac equation in 
the presence a constant magnetic field and a Coulomb potential; we construct a two- 
term mixed-basis in order to compute the energy spectrum. In Sec. 3, we solve the 
Pauli equation and show how the energy spectrum depends on the spin orientation. 
In Sec. 4, with the help of a mesh method we solve the Dirac equation numerically. 
We compare the numerical result with those obtained in the non relativistic limit 
and with the the help of the variational method. Finally we briefly discuss the 
validity of our approach, and further applications. 
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2. Dirac equation 

The covariant Dirac equation 7,20 in atomic units, h = M = e = 1, in the presence 
of an external electromagnetic field can be written as follows 

7^-r,-^) + c)* = o, (i) 

where the Dirac matrices 7 M satisfy the commutation relation {7^, 7"} + = 2g M ' y .Thc 
metric tensor g a $ written in polar coordinates (t, p, i?) has the form: 

g a p = diag(-l, 1,/r 2 ), (2) 

and r p are the spinor connections 21 . The vector potential A a associated with a 
2D Coulomb potential and a constant magnetic field interaction is: 

^ = (-f,0,-^). (3) 
P z 

Choosing to work in the diagonal tetrad gauge we have that the 7 matrices take 
the form 

l p = 7\ 7* = -, f =7°, (4) 
P 

where 7 M are the standard constant Dirac matrices, satisfying the anticommutation 
relations {7^,7"} = 1rf v . 

In the diagonal tetrad gauge, 20 the spinor connections take the form 

r = 0, r 1 = 0, r 2 = i 7 y. (5) 

Since we are working in 2+1 dimensions, it is convenient to introduce the following 
representation of the Dirac matrices in terms of the Pauli matrices 

7 1 = (7i, 7 2 = cr 2 , 7° = CT 3- (6) 

The Dirac equation (1) expressed in the diagonal tetrad gauge commutes with the 
operators ijj and —i-§$, therefore the spinor ^> can be written as 



®(t, p) = -== exp(-iEt + ik^)^(p), (7) 



with 



where the ^fp factor has been introduced in order to eliminate the spinor contri- 
bution 7^r M in the Dirac equation (1). Substituting Eq.(7) into Eq.(l) and taking 
into account Eq. (5), we obtain 



E Z . 2 , k# Bp , , 

— + — +l(T 2 d p +(7 1 (— + — - +(7 3 C 

c cp p 2c 



V> = 0. (9) 
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After substituting the Pauli matrices into Eq. (9) we reduce our problem to the 
following system of coupled differential equations 

^ + (^ + ^W(- + C +-)*i = o, (10) 

dp \ p 2c / c cp 

d ^-(^ + ^)^-(--c + ^=0. (11) 
dp \ p 2c / c cp 

In order to analyze the nature of the eigenvalues k$ of the — iJ^ operator, we 
notice that the spinor \I/ expressed in the (rotating) diagonal tetrad gauge, is related 
to the Cartesian (fixed) spinor ^ c by means of the transformation 20 : 

* = S^)" 1 ^, (12) 

where the matrix S($) can be written as 

= cxp(-^ 7 1 7 2 ) = cxpH^ 3 ). (13) 
Noticing that S(fl) satisfies the relation 

S(d + 2tt) = -5(?9) (14) 

we obtain 

*(<? + 2tt) = -*(!?), (15) 

and we have k$ = N + 1/2, where N is an integer number. The rotating Dirac 
spinor 4" can be written in terms of *5f c as: 

Taking into account the relation (12) between and ^ c as well as (7) we readily 
obtain 



The system of equations (10)-(11) does not admit exact solution in terms of 
special functions, therefore we proceed to compute energy eigenvalues with the 
help of variational techniques. 

In non relativistic quantum mechanics, variational methods provide a power- 
ful technique for the construction of approximate eigenvalues and eigenfunctions 
and for calculations involving sums over the complete energy spectrum. The vari- 
ational method cannot be straightforwardly extended to the Dirac case because 
the Hamiltonian associated with Eq. (1) is not bounded from below. Any positive- 
energy eigenvalue can collapse into a negative-energy eigenvalue as the basis set is 
increased or as the nonlinear parameters of the basis are varied. 22 

The earliest attempts at the solution of the relativistic basis set problem revealed 
a tendency for any variational approach to a bound-state of positive energy to 
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collapse into the negative-energy region of the spectrum. This effect, which was 
first noticed by Kim, 23 has since become known as variational collapse or finite 
basis set disease. 

Among the possible ways to circumvent the problems related to a finite basis and 
the variational collapse, we have chosen to work with a quadratic Dirac equation: 



dp- 



'E Z 

r | - + - 

cp 



cp z 



k§ Bp 



9 k^a" 1 

C H 7T- 



Ba 3 
~2~7 



i> = (18) 



which can be obtained after squaring Eq. (9). It is worth noticing that Eq. (18) is 
a second order differential equation and therefore it does not present the problems 
associated with the variational collapse. 

In order to solve Eq. (18) we propose a mixed-basis variational approach that 
has been successfully applied in the non relativistic case as well as in computing the 
energy eigenvalues of the (2+1) Klein-Gordon equation. We introduce the following 
trial function 



1p = Clip! + C 2 tp2 = C H lpH + C i>0, 



(19) 



where "01 = il>H and ip2 = ipo are the corresponding hydrogen and oscillator wave 
functions associated with the ground state, co and cjy are constants to be cal- 
culated. It is worth mentioning that our basis is not orthogonal under the inner 
product (ipi | ipj) = J °° tpjtpjdp. Substituting (19) into (18), and performing varia- 
tion on the basis coefficients cj , we readily obtain the following matrix equation: 



dp 2 

f.mB 2 E 2 .. IB 2 2EZ 



Z^ B Z 

+ (—- f4)Aij - —Fij + kvGij + i-Hij 

C AC C 



(20) 



with 



^(r)), G 



ipj(r) 



(21) 



Dij = (A(r) \p 2 \^(r)) , Fa = (^(r) \a%\ ^(r)) . 

Instead of computing all the matrix terms in (21) one can use that ?pH and ipo 
are the corresponding hydrogen and oscillator wave functions associated with the 
ground state 

H\V >= c H E H \^H > +c E \^o > ■ (22) 
Using the expression (22) we reduce Eq. (20) to the following matrix equation: 

QijCj = (23) 
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where the components of the square matrix Qij are 



'■'"-if 



2(E-E H )Z 



C 



li 



¥) 2 - D "(0-^ k ' +G ^ (24 » 



Q 



12 



Eh 

c 



2c 



S + 



2(E - E H )Z 



C 



12 



D 



12 



I) r 



Eo 
c 



2EZ 



-C 2 



(26) 



where S =< ^h\^o > is the overlap between the hydrogen and oscillator wave- 
functions. 

The equation det(Qjj) = gives as a result an algebraic equation for E that 
permits one to obtain the variational estimate for the energy. Among the solutions 
of the secular equation we identify the energy as the lowest positive root. 

The advantage of the two-term mixed-basis variational approach is two-fold: 
First, it gives accurate results in the pure-Coulomb as well as in the strong magnetic 
field regimes and very good estimates for the intermediate region. Second, the 
variational wave function can be calculated in a closed form. We do not have to 
handle large terms variational basis and the problems related to balancing small 
and large components of the Hamiltonian. 

In this paper we are interested in computing the energy spectrum of the ground 
s state. The lowest energy configuration corresponds to a state that in the absence 
of the Coulomb interaction has the electron spin directed opposite to the magnetic 
field. 11 

The solution of the Dirac equation in the presence of a constant magnetic field 
with the electron spin directed opposite to B is 



(27) 



the energy associated with spinor (27) does not depend on the magnetic field 
strength and it is given by the simple expression: 




El 



(28) 



which corresponds to a zero non relativistic energy. 

In order to compute the spinor ipn , solution of the (2+1) Dirac equation in the 
presence of the Coulomb interaction we solve Eqs. (10)-(11) for B = 0. The solution 
can be readily obtained in terms of Laguerre Polynomials that for the s state take 
the simple form 



(2\y+i 
y/T(2u+l) 



P e 



-\ P 



y /1 - E H j& 
^l + E H /c< 



(29) 



I 



I 
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where v and 7 are given by 



The energy of the ground state is 

Eh 




(30) 
(31) 

(32) 



Some coefficients of the matrix Qij are divergent in the hydrogen-oscillator basis. 
In order to circumvent this problem we replace p 1 ' 2 by p x ~ v in Eq. (27). This 
substitution permits one to compute the coefficients Aij and does not sensitively 
affect the computation of the energy energy spectrum. 



3. Non relativistic limit 



In this section we are interested in computing the energy energy spectrum of the 
2+1 Coulomb atom in the presence of a constant magnetic field when one consider 
spin effects. The non relativistic limit of the Dirac equation (1) is given by the Pauli 
equation: 



{P+ _Bxrf 



Z 
P 



D 



(33) 



where s is the spin operator. For a constant magnetic field, we have that s- B can 
be written as 



B = ^a 3 B. 



(34) 



Since we are dealing with a two-dimensional problem, we choose to work in polar 
coordinates (p, The angular operator — ids commutes with the Hamiltonian (33), 
consequently we can introduce the following ansatz for the eigenfunction: 



(35) 



^2w y/p 

Substituting (35) into (33), we readily obtain that the radial function u{p) satisfy 
the second order differential equation 



ld 2 u 
Yip 1 



-(m 2 



1)1 

4 V 



2 „2 



B 2 p 
8c 2 



E 



B 



m- 



B 
''2~c 



(36) 



The solution of equation (36) cannot be expressed in terms of special functions. 11 
Taking into account how the spin couples to the magnetic field in the Pauli equation, 
we can compute the energy eigenvalues of Eq. (36). The energy spectrum of the 
Pauli can be computed in a straightforward manner with the help of the solutions of 
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the Schrodinger equation. The presence of the spin introduces a shift in the energy 
proportional to the magnetic field strength. 

Enonrel = Egchr ± 7T ' (^7) 

2c 

where the upper and lower signs correspond to the a magnetic field oriented along 
or in opposite direction to the electron spin respectively. Since we are interested in 
a configuration where the spin is opposite direction to the spin, which corresponds 
to the lowest energy configuration, we choose the negative sign in Eq. (37). 

The computation of energy spectrum of the (2+1) Schrodinger equation in the 
presence of the a Coulomb potential and a constant magnetic is a problem that 
has been widely discussed in the literature and different analytic and numerical 
techniques are at our disposal. 18 



4. Discussion of the results 

In this section we proceed to compute the energy spectrum of the Dirac equation 
with the help of the methods discussed in Sec. II and III. In order to compare 
the accuracy of the variational approach as well as of the non relativistic limit, 
the numerical computations of the relativistic energy spectra are carried out with 
the help of the Schwartz method 19 , which is a generalization of the mesh point 
technique for numerical approximation of functions. 

This method gives highly accurate results given a thoughtful choice of the ref- 
erence function. For Eqs.(10)-(11) we chose as the interpolation function 

/<w = X><^b <»> 

where 

u{p) = sin[7r(p/^) 1/2 ] (39) 

r m is a zero of u(p), a m is a zero of its derivative, and h is the step of the quadratic 
mesh. The use of this scheme on Eq. (10)-(11) leads to an algebraic eigenvalue 
problem, giving a non-symmetric matrix to be diagonalized in order to obtain the 
energy values. This selection gives as result an overall convergence rate which goes 
approximately ase« iq(- 2 /^) n ) where N is the number of truncations. The accu- 
racy of the technique has been verified by computing the energy spectrum of the 
2D Hydrogen atom 14 ' 15,16 and reproducing the analytic results obtained by Taut 17 
for the excited states. 

Fig. 1 shows the dependence of the non relativistic energy E—c 2 on the magnetic 
field strength B. In order the compare the results with those previously reported 
in the literature, 1 ' 18 we introduce the parameter 7 = B/c. Fig. 1 shows that the 
two-term variational method gives very good results for 7' = 7/(7+1) < 0.7 The 
non relativistic approximation works for all values of the magnetic field strength. 
Table 1 shows clearly that the variational approach is better suited for small and 
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Fig. 1. Energy E — c 2 in atomic Rydberg units of the ls~ state as a function of 7' = 7/(7 + 1). The 
long-dashed line is obtained by numerical methods; the solid line corresponds to the mixed-basis 
variational method. The dotted line is obtained by using the Pauli equation 



intermediate values of 7' but the non relativistic Pauli equation guarantees at least 
four figures for any value of B and gives reliable results in a range of values that 
its not restricted to values B << B cr as the semiclassical approximation. 9 

For weak magnetic helds, Ho and Khalilov 9 have obtained approximate energy 
values with the help of a semiclassical method, their result, written in atomic Ry- 
dberg units, takes the form: 

-1/2 



2 



1 + 



1 



(40) 



c 2 (n+V^-Vc 2 ) 2 . 

Expression (40) gives good energy estimates for 7' < 0.05. In fact, for 7' = 0.05 
one obtains E sem ~ —4.02649, a value slightly higher than that one obtained nu- 
merically and with the help of variational techniques. For strong magnetic fields 
the authors obtain recursion relations that determine the coefficients of the se- 
ries expansion of the wave function, the possible energies, and the magnetic fields, 
showing that the two dimensional Dirac electron in the presence of an homogc- 
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Table 1. Relativistic energy values E — c 2 in atomic Rydberg units for s state and a comparison 
with the non relativistic energy spectrum. The first column corresponds to the value obtained via 
the mixed-basis variational method, the second column is the numerical value relativistic energy 
E — c 2 , the third column corresponds to the energy obtained after solving the Schrodinger, the 
fourth column is the energy spectrum of the Pauli equation, the fifth column (diffl) is the difference 
between the the first and second columns, the sixth column (diff2) is difference between the second 
a fourth column. 



7 variational numerical Schrodinger Pauli diffl difF2 



0.01 


-4.010308476 


-4.010303938 


-3.9999905 


-4.01009151 


-4.5385 10~ 6 


0.000212427 


0.05 


-4.052584016 


-4.052582214 


-3.9997404 


-4.052371979 


-1.80165 10~ 6 


0.000210235 


0.10 


-4.110158292 


-4.110161852 


-3.9988434 


-4.109954511 


3.55978 10~ 6 


0.000207341 


0.15 


-4.17375409 


-4.173760085 


-3.9970851 


-4.173555688 


5.99508 lO" 6 


0.000204397 


0.20 


-4.244337352 


-4.244360501 


-3.9941592 


-4.2441592 


2.31488 10~ 5 


0.000201301 


0.25 


-4.32313341 


-4.323172984 


-3.9896419 


-4.322975233 


3.95743 10" 5 


0.000197751 


0.30 


-4.411646466 


-4.411703859 


-3.9829381 


-4.411509529 


5.73929 10~ 5 


0.00019433 


0.35 


-4.511776714 


-4.511854721 


-3.9732025 


-4.511664038 


7.8007 10~ 5 


0.000190683 


0.40 


-4.625971588 


-4.626065158 


-3.9592116 


-4.625878267 


9.35697 10~ 5 


0.000186891 


0.45 


-4.757387992 


-4.75752472 


-3.9391601 


-4.757341918 


0.000136728 


0.000182802 


0.50 


-4.91019842 


-4.910498134 


-3.9103193 


-4.9103193 


0.000299714 


0.000178834 


0.55 


-5.089955586 


-5.090843552 


-3.8684466 


-5.090668822 


0.000887966 


0.00017473 


0.60 


-5.303980092 


-5.306876432 


-3.806706 


-5.306706 


0.00289634 


0.000170432 


0.65 


-5.561460262 


-5.570890138 


-3.7135808 


-5.570723657 


0.009429876 


0.000166481 


0.70 


-5.871135116 


-5.902019599 


-3.5685234 


-5.901856733 


0.030884483 


0.000162866 



neous magnetic field is a quasiexactly solvable problem. Unfortunately, just like in 
the non-relativistic 17 and Klein-Gordon 15 cases, the quasiexactly solvable energies 
do not correspond to the ground state. 

Since the energy spectrum for the ground state, in the pure magnetic field 
regime does not depend on B (28), the variational energy slightly deviates from 
the numerical value for strong magnetic fields. Regarding the methods and results 
presented in this article we have to emphasize that we do not require any special 
identification between small and large components of the Dirac spinor, a condi- 
tion that is present in most of the alternative variational methods available in the 
literature. The electromagnetic configuration given by Eq.(3) presents the interest- 
ing feature of combining a Coulomb potential with a magnetic field in a way that 
makes unsuitable a single trial-function approach. The accuracy of the mixed-basis 
approach can be improved when one uses a many-terms basis, nevertheless this one 
introduces higher order Qij matrices and further complications in the computation 
of the energy eigenvalues. 

Regarding the range of validity of variational and non relativistic approximation, 
we have that both methods give very good results for values of 7 < 1, which 
correspond, for the hydrogen atom, to B < 2.35 x 10 9 G, a value that is far from 
the values of B that one create in a laboratory 24 . 

However strong fields can be mimicked in nanomaterials like GaAs or InSb with 
very small effective masses, where the methods and techniques of quantum effects 
in the presence of strong fields can be applied 24 . In this physical scenario the 
techniques developed in this article could be of help. 
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